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ABSTRACT 

The composition of the overall spectrum of cosmic rays (CRs) is studied under the assumption that 
ultra high energy CRs above the energy 10^^ eV are produced at the shock created by the expanding 
cocoons around active galactic nuclei (AGNs). It is shown that the expected CR composition is 
characterised by two peaks in the energy dependence of the mean CR atomic number < A[e) >. 
The first one at the energy e « 10^^ eV corresponds to the very end of the Galactic CR component, 
produced in supernova remnants (SNRs). It is followed by a sharp decrease of < ^(e) > within 
the energy interval from 10^^ to 10^^ eV. This is a signature of the transition from Galactic to 
extragalactic CRs. The second peak, with < InA >S3 2, at energy e « 10^^ eV, expected at the 
beginning of the GZK cutoff, is the signature of the CR production by the nonrelativistic cocoon 
shocks. The calculated CR composition is consistent with the existing data. The alternative scenario, 
which suggests reacceleration increasing the energy of CRs produced in SNRs by a factor of 30, is also 
examined. 

Subject headings: acceleration of particles — cosmic rays — galaxies: active — shock waves — super- 
nova remnants 



1. INTRODUCTION 

The overall origin of cosmic rays (CR), in particular at 
ultra high energies, above 10^^ eV, is still an unresolved 
problem in astrophysics (e.g. iNagano fc WatsonI 120001 : 
lStanevl[2f)0l lCronin][20f)5t llnoud 12001 7" Understanding 
this origin requires the determination of the astrophysi- 
cal objects, that are the CR sources, and of the appro- 
priate acceleration processes, that form the CR spectrum 
in these objects. During the last several years consider- 
able progress has been achieved in this field, both exper- 
imentally and theoretically. Recently the sharp steep- 
ening of the CR spect rum above 3 x 10^^ eV was es- 
tabhshed in the H i Res (Bergman et al.ll2007t l and Auger 
(|Yamamoto et all 12007 ) experiments. It presumably 
corresponds to the so-called Greizen-Zatsepin-Kuzmin 
(GZK) cutoff, caused by energy losses of the CRs in 
their interactions with the microwave background radi- 
ation. This is evidence that the highest energy part of 
the CR spectrum is of extragalactic origin. It was also 
recently demonstrated (|Berezhko fc Volkll2007f ) that the 
CRs with energies up to e ~ 10^^ eV can be produced in 
supernova re mnants (SNRs) as a result o f diffusive shock 
acceleration (jKrvmskvl [l977t iBelll I1978D . and that the 
observed CR energy spectrum can be well represented 
by two components: the first, dominant up to 10^^ eV, 
consists of CRs produced in Galactic SNRs; the second, 
dominant at energies above 10^* eV, is of extragalactic 
origin. The latter component can be produced at the 
shock that is formed by the expanding cocoon around 
active galactic nuclei (AGNs): above the energy 10^^ eV 
the overall energy spectrum of CRs, produced during 
the AGN evolution and released into intergalactic space, 
has an appropriate power law form, which extends at 
least up to the energy emax ~ 10^° eV, well above the 
GZK cutoff (lBerezhkoll200l . This is the "dip scenario" 



of th e overall CR spectrum formation (|Berezinskv et al] 
|2006[) . Within the alternative "ankle scenario" , extra- 
galactic CRs are expected to dominate only above the 
energy 10^^ eV (Bcrezinskv ct al. 20 ^^ . It was noted 
(jBerezhko fc V61k,2007; .Berezhkg, 20081) . that the com- 
position of CRs at e > 10^'' eV is expected to be very 
different in these two scenarios. Therefore the study of 
CR composition can in principle provide the signature of 
the acceleration mechanism which creates the spectrum 
of extragalactic CR component. 

Here we present the calculations of the expected com- 
position of CRs, accelerated by cocoon shocks and 
demonstrate that it has well pronounced peculiarities, 
which can be considered as a signature of ultra high en- 
ergy CR production by nonrelativistic shocks. The cal- 
culated composition is compared with the composition 
expected within the alternative scenario and with the 
existing data. The alternative scenario, which suggests 
reacceleration increasing the energy of CRs produced in 
SNRs by a factor of 30, is also examined. 

2. RESULTS AND DISCUSSION 

The composition of CRs is determined if we know not 
only their all-particle spectrum. 
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that is the differential intensity J(e) with respect to par- 
ticle energy e, but also the spectra Ja{^) of all rele- 
vant elements with atomic number A. Since the shock- 
accelerated CRs come originally from the thermal back- 
ground plasma, one has to expect that the flux of each 
CR element in/near the source, is proportional to the 
number density of this element A'a in the background 
plasma: J%{e) oc Na- Therefore it is useful to represent 
the spectrum of each CR element in the form 

Ji(e) = eA(e)aAJ|i(e), (2) 

where a a = A^a/A'h is the abundance of element A rel- 
ative to the hydrogen abundance and e^ie) is an en- 
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richment factor which describes the preferential produc- 
tion of element A relative to the production of protons 
(marked by the subscript H). 

Direct measurements of CR fluxes at energies e < 
10^'^ eV show that the factor ca > 1 considerably ex- 
ceeds unity and t hat it progressively increa.ses w ith the 
increase of A fe.g. iBerezhko fc Ksenofontovlll99"9f) . This 
means a considerable enrichment of CRs in heavy ele- 
ments compared with the interstellar abundance. Note 
that at these relatively low energies, where direct deter- 
mination of CR particle parameters is possible, the fluxes 
Ja and Jh are usually compared at the same energy per 
nucleon, that is for ea = Ae^. Since it is impossible to 
determine the mass of individual CR particles at ultra 
high energies, the particle energy instead of the energy 
per nucleon is the most appropriate variable for the anal- 
ysis of CR spectrum and composition in this case. 

The observed CR composition at e < 10^^ eV 
is well consistent with the scenario in which Galac- 
tic CR component is p roduced by supernova shocks 
(|Berezhko fc Volkl |2007|) . The enrichment of shock- 
accelerated CRs by heavy elements is due to a more effec- 
tive injection and subsequent accel eration of heavy ions 
(see IBerezhko fc Ksenofont"ovlll999l for the details). 

The spectra of ultrahigh energy CRs produced by non- 
relativistic coco on shocks can be represented in the form 
(|Berezhkoll20"08h 



Jl(e) = CAe- 



(3) 



where 7 « 2.6 and the values of parameters Ca oc baNa 
are determined by the background medium's composition 
Na and by the enrichment factors ca- Since the injec- 
tion/acceleration of CRs at cocoon and supernova shocks 
are expected to be very similar, the values of the enrich- 
ment factors BA in the former case can be extracted from 
the experimentally measured fluxes of the Galactic CR 
con iponent JA(e) at some fix ed energy, say e = 1 TeV (see 
e.g. IBerezhko fc Volkl[2007[ ). according to the expression 

CA = Ja(1 TeV)/JH(l TeV)(iVH/iVA)0. (4) 

Here {Na/Nb_)q is the solar system relative abundance 
of elements with atomic number A. 

In order to calculate the spectra of extragalactic CRs 
inside the Galaxy one needs to take into account the 
modification of the CR spectra due to their interaction 
with the cosmic microwave background (CMB) and due 
to their modulation in the Galactic wind. This can be 
done by representing the observed extragalactic CR spec- 
trum in the form 



(5) 



where factors Pi and Ph describe CR flux modification 
due to Galactic wind and CMB respectively. Since the 
diffusive mobility of the CRs is inversely proportional 
to the particle charge number Z, the first factor can be 
represented in the form 

Pi = exp(-Zei/e), (6) 

where ei - lO" eV (jBerezhko fc Volkl l2007f ). This ex- 
pression shows that at a given energy the lighter CR 
species can more easily penetrate the inner Galaxy from 
outside. Due to this factor the observed extragalactic 
CR component becomes progressively heavier with in- 
creasing energy. 
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Fig. 1. — Overall CR intensity as a function of energy (thick solid 
line), Galactic CR component produced in SNRs (thick dashed 
line) and extragalactic component produced in the IGM (thick 
dash-dotted line). The spectra of different extragalactic compo- 
nents are shown by thin linos. Experimental data obtained in the 
ATIC-2 (Panov ct al. 2006), JACEE( Asakimori ct al. 2003), KAS- 
CADE tAntoni et aU2005.) . Akeno (.Takeda et al..200a) . and HiRes 
IIBeremarr'er'ani2007l) experiments are shown as well. 

To describe the effects of CRs interacting with the 
CMB radiation during their propagation from extra- 
galactic sources to the observer - pion and pair produc- 
tion, CR nuclei photodesintegration and adiabatic cool- 
ing - we use a simple analytical form of the factor Ph 



Ph = exp[-(^e/eh)']. 



(7) 



At Eh — 1-3 X 10^^ eV for the case of helium and iron 
nuclei it satisfactorily fits the results of detailed cal- 
culations which c onsistently describes all these effects 
(|Berezinskv et al][20 06). In the case of protons, instead 
of the modification factor P h we directly use their m od- 
ified spectrum calculated bv iBerezinskv et all ()2006l ). 

As is clear from expression (7), the interaction of CR 
nuclei with the CMB radiation produces a sharp cutoff 
in their spectra at the energy, which is proportional to 
the nuclear mass. Due to this fact CRs are expected to 
become progressively heavier at energies above 10^^ eV. 

The cocoon shock propagates initially through the in- 
terstellar medium (ISM) of the host galaxy and later 
through the intergalactic medium (IGM). The IGM 
is less abundant in heavy elements than the ISM: 
the median metallicity of the IGM is (Na/N-h) — 
0.2{N\/N yi)(t, for al l relevant elements heavier than he- 
lium (e.g. ICen fc Os triker 2006). Since the maximal 
CR energy produced by the expand ing cocoon shock 
decreases with time (jBerezhko! l2008l ) . the local metal- 
licity at the shock front changes from the ISM value 
{Na/Nh) w (TVa/A'^h)© at early phases, when the shock 
produces CRs up to e^ax ~ 2^ x 10^" eV, to the IGM 
value at late phases, when Cmax ^ Z x 10^* eV. To 
describe this effect we uses a metallicity (Na/Nh) = 
O.2[e/(ZxlOi8eV)]°-35(7VA/iVH)0. For helium we use the 
cosmological value A^hc = O.OSTVh (e.g. IBerezinskv et al.l 
200 6). Using these number densities one can calculate 
the relative values of the coefficients Ca/Ch for all el- 
ements except protons. The proton coefficient Ch de- 
pends on a number of physical factors such as the IGM 
hydrogen number density A^h and the power and spatial 
distribution of nearby AGNs. Its value is determined by 
the fit of the expected all-particle CR flux J(e) to the 
existing measurements of the CR intensity at energies 
e > IOI8 eV. 



Composition of CRs from AGNs 



3 



'I I I 



.ATIC-2 ♦JACEE 




Fig. 2. — Mean logarithm of the CR nucleus atomic number as a 
function of energy. Solid and dashed lines correspond to the dip and 
ankle scenarios respectively. Experimental data obtain ed in the 
ATIC -2. JACEE, KASCADE (QGSJET and SYBYLL: IHorandell 
[200ah . and HiRes (QGSJET and SYBYLL; lAbbasi et al. 1 120051) 
experiments are shown. 

Secondary protons appearing as a result of photodesin- 
tegration of helium nuclei are also included in an approx- 
imate way: it is assumed that the secondary protons 
(whose total number equals four times the number of 
helium nuclei in the source spectrum with energy above 
2 X 10^* eV) have the normal energy distribution around 
eo = 5 X 10^^ eV within the interval Ae = cq. 

In Fig.l we present the all-particle spectrum which 
includes two coinpone nts: CRs produced in SNRs 
(jBerezhko fc V6l3 l2007t ) and extragalactic CRs, which 
consist of protons (H), helium (He) and three groups of 
heavier nuclei, produced in AGNs. Note, that the dip at 
energy e « 10^^ eV is still clearly seen even though it is 
less pronounced compared with the case of pure proton 
composition. The calculated overall CR spectrum is in a 
satisfactory way consistent with the experimental data, 
except in the energy interval 10^^ < e < 10^^ eV in the 
transition region, where, contrary to the experiment, the 
calculated spectrum has a small dip. Taking into account 
the existing experimental uncertainties at these energies 
it is not clear whether this discrepancy is a serious prob- 
lem for the model or not. Note, that the size and the 
shape of this peculiarity is very sensitive to the details 
of the galactic CR high energy cutoff and to the extra- 
galactic low energy spectrum part. Therefore it is hard 
to predict it quantitatively. At the same time the overall 
spectrum, which is a mixture of two independent compo- 
nents, should have such peculiarity within the transition 
region. If the actual CR spectrum is indeed so uniform 
within the energy range 10^^ < e < IQis eV, as it looks 
according to the Akeno data, then it should be consid- 
ered as a serious problem for the dip scenario. 

The mean logarithm of CR atomic number is repre- 
sented in Fig. 2 as a function of ene rgy. As already quali- 
tatively predicted (jBerezhkdl200l ). the CR mass energy 
dependence < \nA{e) > has two peaks. The first one 
at the energy e « 10^^ eV corr esponds to the very en d 
of the galactic CR component (jBerezhko fc V6l5l2007t ). 
whereas the second, at the energy e « 10^^ eV, is at the 
beginning of the GZK cutoff. Note that at high ener- 
gies e > 10^^ eV information about CR composition is 
obtained from the mean values of the shower maxima 
^max (usually measured in g/cm^), determined by the 
ground based detectors. Knowing the average depth of 



the shower maximum for protons, X^^^, and for iron nu- 
clei, A^j^axJ from simulations, the mean logarithmic mass 
can be estimated from the measured X^ax according to 
the relation (e.g. Horandcl 2005) 

< In A >= (X,„ax - X^,,J/{Xll, - XH^J < In 56 > . 

This conversion requires the choice of a particle interac- 
tion model. Here we use the values of Xmaxi measured 
in the HiRes experiment (jAbbasi et al. l[2005f ). and the 
QGSJET and SYBYLL models to determine <lnA>. 
It is seen from Fig. 2 that the HiRes data are very well 
consistent with the expected sharp decrease of < In A > 
within the energy interval lO^-lO^^ eV. At higher ener- 
gies, between 10"'^* and 10^^ eV, the experimental values 
of < \nA{e) > have a quite irregular behavior. Nev- 
ertheless, the experiment reveals a trend of progressive 
increase of the mean CR mass so that the expected peak 
value < In A >« 1.7 at the energy e ^ 10^^ eV is consis- 
tent with the existing HiRes data. It is also clear from 
Fig. 2, that a more precise experimental determination of 
the CR mass composition for e > 10^^ eV is required for 
a more complete conclusion about the CR origin. 

3. ALTERNATIVE SCENARIO 

The alternative scena rio for the overall CR spectrum 
is the "ankle scenario" (|Berezinskv et al]|2006l ). In this 
case the extragalactic source spectrum, as compared 
with the dip-scenario, is assumed to be much harder 
oc so that it becomes dominant above an en- 
ergy of e = 10^^ eV (e.g. iBerezinskv et al.ll2006l) . There- 
fore, to fit the observed overall CR spectrum one needs a 
third component to fill the gap between the Galactic CR 
spectrum, produced by SNRs, and the hard extragalac- 
tic spectrum. It is hard to believe that this third CR 
component could be completely independent of the CR 
component produced in SNRs, since in such a case one 
should expect some peculiarity in the spectral shape at 
energy e « lO^'' eV, where these two components are ex- 
pected to match. The appropriate solution of this prob- 
lem is the existence of some kind of reacceleration process 
which picks up the most energetic CRs from SNRs and 
substantially increases their energy, resulting in a smooth 
extension of the first CR component towards the higher 
energies. 

We model the spectra of reaccelerated CRs in the fol- 
lowing way, without specifying the reacceleration mecha- 
nism. For every element with the nuclear charge number 
Z, as in the dip scenario, we use the spectrum, which 
coincides with Ja{^) for e < Cmaxi ^ form 

JA(e) = JA(eLxi)(e/emaxi)~^exp(-e/e^ax2) (9) 

at e > e^axi ■ Here e^^xi is the minimum energy of parti- 
cles involved in reacceleration, and e^ax2 the maximum 
particle energy achieved during reacceleration. It is nat- 
ural to assume that these energies scale proportional to 
the rigidity e^^x — ^^max- Here the superscript p de- 
notes protons. The quantities Cmaxii ^max2 ^^'^ 7 
treated as free parameters whose values are determined 
as a result of the best fit. CR acceleration by spiral den- 
sity shocks in the Galactic Wind (Volk & Zirakashvi3 
2004 ) , or acceleration in the vicinity of pulsars (Bell] 
1992t lBerezhko1ll994[ ) could play the role of the reaccel- 
eration mechanism. 
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Fig. 3. — The same as in Fig.l, but for the ankle scenario. 
The dashed hne represents the Galactic component, which includes 
CRs produced in SNRs and reaccelerated CRs. The dash-dotted 
line represents the extragal actic component; it cor responds to the 
source spectrum J\ oc l)Berezinskv et aLll2006I V 

We present in Fig. 3 the CR spectrum calculated within 
the ankle scenario, with the adopted values e^^xi ~ 
5 X 10^^ eV, eP^^2 = 1-5 x lO" eV and 7 = 3, which 
provid e the best fit to the obs erved CR spectrum. Fol- 
lowing |Berezinsk^^r^n (|2006l ) we use the extragalactic 
CR source spectrum in the form J\ oc and assume 
the composition of these CRs with < ^(e) >= 1. 

It is seen that at e > 10^^ eV it is well consistent with 
the data. However, such a well-known peculiarity in CR 
spectrum as the knee at e « 3 x 10^^ eV is much less 
pronounced in the theoretical CR spectrum than in the 
experiment. This may be considered as an indication 
against the ankle scenario. 

As it is seen in Fig. 2, the CR composition correspond- 
ing to the ankle scenario is considerably heavier at en- 
ergies lO-'^^e < 10^^ eV than in the dip scenario and is 
inconsistent with the HiRes measurements. Note, that 
CR composition of the reaccelerated CRs is a direct con- 
sequence of the heavy CR composition at the energy 
e > e^axii where the reacceleration is assumed to start, 
and it is not sensitive to the details of reacceleration pro- 
cess. 

Since the expected spectrum of CRs produced in 
gamma-ray bursts (GRBs) has a form J" oc with 
7 — 2 — 2.2 only ankle scenario is possible if GRBs 
are the the source of ex tragalactic CR component 
(see iBerezinsky et al.l I2006L and references therein) . It 



means that HiRes CR composition is against GRBs as 
a source of ultr a high energy CRs. In addition, as 
iBerezinskv et al.l ()2006[) argued, the energy output of 
GRBs has a serious problem if they are considered as 
the main source of the extragalactic CRs. Note also, 
that based on the data collected at the Auger experi- 
ment, correlation between the arrival directions of CRs 
with energy above 6 x 10^^ eV and the pos ition of nearby 
AGNs has been find (jAbraham et al.ll2007l ). that strongly 
supports AGNs as a prime candidate for the source of ul- 
tra high energy CRs. 

4. SUMMARY 

The composition of ultra high energy CRs produced 
by nonrelativistic cocoon shocks around AGNs is char- 
acterised by well pronounced peculiarities which are two 
peaks in the energy dependence of the mean CR atomic 
number < A{t) >. The first peak at the energy e ~ 
10^^ eV corresponds to the ver y end of the Galactic CR 
component, produced in SNRs (jBerezhko fc Vol3l2007f l. 
whereas the second, at the energy e w 10"^^ eV, is ex- 
pected at the beginning of the GZK cutoff. The strong 
energy dependence of the CR composition within the 
energy interval from 10^^ to 10^^ eV is expected as a 
signature of the transition from heavy Galactic to light 
extragalactic CRs, whereas the detection of a heavy CR 
composition at e « 10^^ eV has to be considered as the 
signature of CR production by the nonrelativistic cocoon 
shocks. 

The existing measurements of CR composition arc con- 
sistent with the dip scenario with a formation of the CR 
spectrum in SNRs and AGNs and are inconsistent with 
the ankle scenario, which includes reacceleration of CR 
produced in SNRs as a third CR component. Additional 
peculiarity in the overall CR spectrum - the dip in the 
transition region 10^^ < e < 10^^ eV - is expected within 
the dip scenario. Since it is not seen in the existing data 
it is a real difficulty for the dip scenario. It is therefore 
clear that a more precise measurements of CR spectrum 
and composition at energies above 10^^ eV are needed 
for a strict determination of CR origin. 
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